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Abstract. In this paper we consider vacuum Kasner spacetimes, focusing on 
those that can be parametrized as linear perturbations of the special Petrov type 
D case. For these quasi-D Kasner models we first investigate the modification 
to the principal null directions, then a Teukolsky Master Equation for fields of 
any spin, considering in particular the quasi-D models as curvature perturbations 
of the type D background. Considering the speciality index and the principal 
null directions and comparing the results for the exact solutions and those for the 
perturbative ones, this simple Kasner example allows us to clarify that perturbed 
spacetime do not retain in general the speciality character of the background. 
There are four distinct principal null directions, although they are not necessarily 
first order perturbations of the background principal null directions, as our 
example of the quasi-D Kasner models shows. For the quasi-D Kasner models 
the use of a Teukolsky Master Equation, a classical tool for studying black hole 
perturbations, allows us to show, from a completely new point of view, the well 
known absence of gravitational waves in Kasner spacetimes. This result, used 
together with an explicit expression of the electric Weyl tensor in terms of Weyl 
scalars, provides an example of the fact that the presence of transverse curvature 
terms does not necessarily imply the presence of gravitational waves. 



1. Introduction 

In this paper we apply, to the vacuum Kasner spacetime models [1] , some mathematical 
tools developed in the context of black hole physics, i.e. the speciality index (SI) [2] 
and the Teukolsky Master Equation (TME) [3]. The first one indicates whether 
a given spacetime is algebraically special and is used here, in conjunction with a 
study of the principal null directions (PNDs), to discuss how linear perturbations 
may change the speciality condition. Kasner models are Petrov type / (general), yet 
they are simple enough that a direct comparison between the exact and perturbative 
solutions is possible, which is very useful to clarify, in a transparent way, some results 
of perturbation theory. In particular, the usually defined SI S only changes at second 
order for perturbations of a generic type D background [2], yet one would expect 
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the perturbed spacetime to be of general type I Petrov type. Considering quasi-f 
Kasner models that can be parametrized as linear perturbations of the special Petrov 
type D case, we show that they have four distinct PND, as expected and in complete 
agreement with the exact models. However, the PND of a linearly perturbed model 
arc not necessarily first order perturbations of the background PND, as we show for 
our quasi- Kasner models. 

The TME is a single partial differential equation describing the behavior of 
perturbations of a given Petrov type D background due to fields of any spin 
(0, ±1/2, ±1, ±3/2, ±2); it is therefore a classical tool for studying perturbations of 
black holes, in particular gravitational radiation [3]. It is used here to show the 
absence of gravitational waves in the quasi- 1? Kasner models. This confirms, with a 
completely new approach, the well known absence of gravitational radiation in general 
type I Kasner spacetimes, shown for example by the vanishing in Kasner models of 
the gravitational superenergy flux (see e.g. [4]). Then, we use the non vanishing 
of the Weyl scalar ipQ and ipi in a general Petrov type / Kasner model (either an 
exact or pertubative solution) to illustrate that the non- vanishing of these transverse 
contributions to the electric Weyl tensor (the only curvature contribution in vacuum 
to the Jacobi or geodesic deviation equation) does not necessarily imply the presence 
of gravitation radiation. 

As the Ncwman-Penrosc formalism will be adopted here, conventions and notation 
will follow the existing literature on the subject [3,5], including that for the metric 
signature (+,-,-,-). 

2. Petrov classification: speciality index and principal null directions 

Defining the complex tensor Cabcd = Cabcd — i*Cabcd, one can introduce the two 
complex curvature invariants, in tensor and Newman-Penrose (NP) form 

/ = ^CabcdC"'''' = (V'0V4 - 4VlV3 + 3V^) (2.1) 

and 

These can be used to define the speciality index [2,6] 
27 /2 

5 = ^; (2.3) 

with its value this demarcates, in an invariant way, the transition from certain 
algebraically special solutions (5 = 1) and the general Petrov type / {S ^ 1) [7]. 
This quantity can be easily evaluated in the case of the vacuum Kasner [8] metric: 

ds^ = dt^ - t^P^dx"" - t^P'dy^ - e^^dz^ , (2.4) 

where 

Pi+P2+Ps=pI+pI+pI = 1 (2.5) 
and for the moment there is no preferential ordering of the Kasner indices. 
Let us now introduce a Newman-Penrose tetrad 

/ =l^[d,+t-P'd,], n=^[dt-t-^^d,\, (2.6) 
m=^[t-^^dy + it-^-d,], (2.7) 
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which gives the non zero spin coefficients: 

,= ^^ip.+P,)^^P, ^-^ = -7, A=^(p2-P3) = -a, (2.8) 
the non zero Weyl scalars: 

, , Pl{P2-P3) , P2P3 

V'o = V'4 = ^ , "^2 = --^, (2.9) 

and generates the time independent SI 
27 

S = -plil-p,). (2.10) 

This tetrad is usuaUy caUed "transverse" because of the property tpi — 'ipz — [6,9], as 
it is expHcit from equation 12.14|l below (see also the discussion at the end of Sections 
3 and 4). From H2.10|l clearly the Petrov type is / in general; the null tetrad may be 
said in canonical form because tpo = ip4 [10]. 

From H2.10|l S results well defined for any Kasner solution, while in general this 
may not be the case [6] . The Kasner metric admits two special subcases when two of 
the Pi indices are equal: it then follows from (|2.5(l that either pi = p2 — 0, p3 — 1 
(and permutations) and the spacetime is flat in this case, or pi = —1/3, P2 — Ps = 2/3 
(and permutations) and one has the Kasner LRS type D solution, with 5=1, with a 
spindle-hke singularity [1,7]. As it was shown by Geroch [11,12], this type D model 
can be obtained as the limit for A/ ~> cx3 of the Schwarzschild solution, provided a 
specific coordinate transformation. For the type D case with p2 — Ps = 2/3 the null 
tetrad above is a principal one, with I and n aligned along the PNDs of the Weyl 
tensor in this case, and is also canonical [10], with ■00 = V'4 = 0- The same tetrad is 
not canonical for the other two physically equivalent type D cases pi = p2 = 2/3 and 
Pi = P3 = 2/3, with ipo = ip4 ^ 0. 

The orthonormal frame naturally associated to H2.6|l - H2.7(l is 

eo = -^{1 + n) ^dt^u, ei = -^{l - n) = t^P^d^; , 
e2 = ^im + fh) = t-P^dy, 63 = ^{m~fn)=t-P^d^. (2.11) 



V2" ' ' V2t 

It is adapted to the static preferential observers with 4-velocity u who use the 
Killing vectors dx,dy,dz to build their spatial axes and therefore directly observe 
the homogeneity of the spacetime. They also observe a purely electric Weyl tensor 

E{u) = —j^ ei (g) ei + 62 (g) 62 + 63 (g) 63, (2.12) 

H{u) = 0, (2.13) 
the electric and magnetic part of the Weyl tensor being 

respectively. It is also useful to give here the expression of E{u) in terms of NP 
quantities in a transverse NP frame 

E{u) = Re(V'2)6c - ^Re(V'o + i'i)eT+ + ^Im(V'o - V'4)6tx , (2.14) 

where ec = 61 (g 61 + 62 (8)62 — 263(8)63, exx — 61(862 + 62(861 and 6t+ = 61(861 — 62(862 
respectively represent a Coulombian and two transverse basis tensors. This relation 
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is further simplified if the transverse frame is also canonical (as the frame H2.6|l - H2.7|l 
used in this case) , giving the result (|2.12|) . 

We now want to compare the above results for the Petrov classification of 
the Kasner exact solutions and results for Kasner models in a neighborhood e [e 
being a small quantity) of the type D case, using p^ = 2/3 + e. Guided by 
the results for the exact solutions, we will be able to interpret the perturbative 
results in a coherent framework. Expanding the SI around the type D case, one 
gets 5 = 1 — 27/4e^ — 27/4£^, which confirms the perturbative result found in [2] 
for a generic type D background starting from the canonical tetrad for type D 
{ip^^ = ■tp^'' = ip-}^^ = V'i"'' = 0)i with the index N = 0,1,... in -tp^^ denoting 
the perturbative order: 

5 = l-3^^e^ + o(e3). (2.15) 

From this result, solutions like this type D Kasner one could seem to change 
perturbatively their Petrov type at second order only. Moreover if the gravitational 
perturbations are algebraically special [5,13] {i.e. either ip^^'' or ip'^'^ vanishing) the 
change should start at third order. One is therefore tempted to conclude from (|2.15|l 
that a spacetime that is a linear perturbation of a type D background retains the type 
D character at first order [14]. On the other hand, it seems very difficult to believe 
that linear perturbations of an algebraically special solution do not alter the speciality 
conditions, i.e. do not modify the number of the PNDs. In other words, the fact that 
S only changes at second order seems to have more to do with its definition that with 
a real speciality of the Petrov type of the linearly perturbed spacetime. Indeed, one 
can redefine the speciality index using a monotonic function of S in order to have 
perturbative changes starting at first order in e: this can be easily shown by using for 
instance the normalized index iSnorm = — 1|- Thus the point is, that once first 
order dynamical variables such as the perturbed metric are chosen, derived quantities 
such as S and iSnorm may or may not be linear in the perturbations, a fact that only 
reflects the way they are defined. 

In order to completely clarify this issue of the Petrov classification of perturbed 
solutions, we are now going to build up, as a concrete example, the principal 
null directions of the general Kasner vacuum models (|2.4|l and to expand them 
perturbatively around the type D case P2 = Ps = 2/3. Of course, as soon as p2 — Ps 
is nonzero (no matter how small) it is clear from H2.9|l and (|2.10|) that the Petrov 
type of any Kasner solution is general, so that there exist four distinct PNDs. As a 
consequence, we will conclude in particular that there must be four perturbed distinct 
PNDs for the quasi- 1? Kasner model, as then confirmed by our perturbative analysis. 
We proceed as follows. First of all, let us recall that, in general, once one has found 
a generic NP frame with its associated Weyl scalars, given in our case by equations 
(|231)-(E7|) and 1121), the PNDs of that manifold are [5] 

hf^^ oc r + bl^^m^ + b(j)m'' + bf^j^bl^^n" , (2.16) 
where the are the four roots {j = 1, 2, 3, 4) of the equation (assuming -04 7^ 0) 

ipo + 46?Ai + + 46^3 + &V4 = . (2.17) 
In our case they result in 

EE = C^Jc + rjVc^ , C = ^^'P' , (2.18) 

" Pl(P2-P3) 
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where C and rj assume independently the values ±1. Assuming, with no loss of 
generality, the following ordered parametrization for Pi,P2,P3 [7,8] 

-^<Pi<0, 0<P2<| ^<P3<1, ^C>3, (2.19) 

the four roots H2.18|l result real, simplifying considerably the calculations. We stress 
that the PNDs H2.16|l have still the freedom to be rescaled by a multiplicative factor, 
which we will fix in order to recover the correct / and n in (|2.6|) in the type D limit. To 
this purpose let us study first the case rj = 1 (still maintaining the freedom ^ = ±1) 
in (|2.18|l . denoting it as In this case, the first two PNDs, using (|2.1t)|) and 

lE!Sll-(EI3), result in 

1 



"(+) = l2— 
"(+) 



(2.20) 



For pi = —1/3, P2 = P3 — 2/3 these directions give exactly the type D limit of n in 
equation H2.6|l . which is a double principal null direction for the spacetime. Let us 
study now the case rj = —1 in (|2.18() . denoting it as In this case the remaining 

two PNDs, using ifTTH)! and ^^-^J^, are 

= r + b^^^m" + 6(_)m'' + (2.21) 

where we point out the absence of the rescaling with respect to the previous PNDs. 
In the same limit as above, we obtain exactly the other repeated PND I of the type D 
case in (|2.6ll . It can be easily verified that in general these four PNDs are not linearly 
independent because the z component is always zero. This is expected because of the 
vanishing of the magnetic part of the Weyl tensor (see ref . [7] , page 54) . 

Note that although the roots in H2.18|l blow up in the type D limit, while 
the 6(_) vanish, the PNDs above are always well defined. However, although the 
type D limit of all the PNDs exists, they don't admit a Taylor expansion around 
this point. With no loss of generality, we demonstrate this property in the simpler 
case of the ordered range (|2.19ll . Assume again = 2/3 + e, with e > to satisfy 
H2.19|l . The perturbation of Kasner constraints (|2.5|1 leads to pi = — 1/3 + O(e^) and 
P2 = 2/3 — e + Ole^). Although the null vectors l,n,m of H2.6|l - H2.7|l have smooth 
Taylor expansions in e, for the PNDs H2.20|l and H2.21|l we obtain 

^ \l + 0{e),ti +0{e), -Ct-h^ +0{ei),o\ , (2.22) 



-1- 



l + 0(e), -t3 +0(e), -Ct-3e5 +0((r5), , (2.23) 



clearly not a Taylor series, but nonetheless manifesting the splitting of the PNDs in 
passing from the Petrov type D to the type /. This peculiar phenomenon regards the 
y component of the null directions, while the z component of these PNDs remains 
always zero. Analogous results will be obtained using another choice of order for the 
Kasner parameters. The nonanalytic behaviour of the principal null directions in the 
parameter space is to be expected at the algebraically special loci. Small changes in 
the parameters have the effect that the coincident pairs of PND bifurcate. We remark 
that if we expand a generic type / Kasner around any point different from the type 
D one instead a Taylor series always exists. 
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3. Teukolsky Master Equation 

Since the Kasner vacuum spacetime with indices pi = — 1/3, p2 = P3 = 2/3 is of Petrov 
type D, we can export the discussion typical of perturbed black hole spacetimes into 
the arena of cosmology, applying the Teukolsky Master equations machinery [3] . 

To study curvature pertubations here it is convenient to introduce another 
Newman-Penrose tetrad, obtained from (|2.6I) - (|2.7|I after a class III null rotation (a 
boost) which makes e vanishing. This is 

I ^i^[dt + t-p'd,] , 71 = [tp'dt - d,] , 
1 

72 

which gives the non zero spin coefficients 

1 1 1 



m = — [t-P'dy+it-P^d,] , (3.1) 



(i-Pi)iPi-i = -<2pv, j = -^pitP'-' 



^ = -^iP^-P^)t'''' = '^^ (3-2) 
and the non zero Weyl scalars: 

^-2(pi+l) I 
V'O = ^ Pl{P2 - Ps) = t~^P^1pA , -02 = - 2*~^P2P3 ■ (3.3) 

The advantage of this tetrad is that it results again principal and canonical [tpQ — 
ip4 = 0) for the Petrov type D with p2 — Ps — 2/3; moreover it satisfies the useful 
condition on the spin coefficients e = (it is a Kinnersley frame) which allows us to 
apply in toto the original Teukolsky notation. In particular the fact that the tetrad 
above becomes the canonical one for the Kasner type D solution with p2 = Ps = 2/3 
ensures the gauge and tetrad invariance of the perturbed ipQ and ip^ [3, 15] for the 
quasi-Z? models. More in general, any covariantly defined quantity that vanishes in 
the background spacetime is a gauge invariant perturbation [15-17]. 

The equations for the gauge and tetrad invariant first order massless perturbations 
of various spin-weight s (coincident with the helicity [18]) in this background, written 
for simplicity in absence of perturbative sources using the frame (|3.1|) and its derived 
NP quantities H3.2|) and H3.3(l . are given by the following NP relations valid for any 
vacuum type D geometry 

{[D- p* +e* +6- 2s{p + e)]{A + fi- 2s-f) (3.4) 

-[5 + n* - a* +P- 2s{t + f3)] {5* + tt - 2sa) - 2{s - l)(s - 1/2)V'2}V'^''^ = 

for spin weights s = 1/2, 1, 2 and 

{[A - 7* + - 7 - 2s(7 + n)]{D-p~ 2se) (3.5) 

-[S* -T* +f3* -a- 2s{a + 7r)](5 - t - 2s/?) - 2(s + l)(s + l/2)7/;2}V''^^ = 

for s = —1/2, —1, —2. The case s = ±3/2 is available in the literature in Geroch-Held- 
Penrose form [19,20] and finally the case s = is given by [21] 

[DA + AZ) - 5*5 - 55* + (-7 - 7* + + A^*)^ + {e + e* - p* - p)A 

+ {~fr - TT + a + T*)5 + (-TT* +T- 13 + a*)5*]'ip^^'^ = . (3.6) 
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Perturbative sources can be easily added following Teukolsky. Introducing a 
"connection vector" that for our Kasner models has the components 

r* = \t-^ , = -t-i , = = , (3.7) 

such that V^Fa — and FTa — 4-02, all these equations collapse to the unique PDE 
form [22-24]: 

[(V- + sr-){Va + sTa) - 45^^]^^^^ - (3.8) 

where ip^"^ is the background Weyl scalar in (|3.3(l . Equation (|3.8(l gives a common 
structure for these massless fields in the LRS Kasner background varying the spin index 
s. The components of the TME with negative spin must be multiplied by a certain 
prefactor in order to give the physical components of the fields, as well explained 
by Teukolsky [3]. For the aim of our subsequent analysis we shall be interested in 
the gravitational case only: consequently one has to consider the solution ip^^") of the 
TME in the case s — ~2; the Weyl scalar with negative helicity is then given by 
1^4 = p'^'il)^~'^\ We recall that this approach is possible for the special type D (and 
type O) spacetimes only, because it is this speciality type that allows to decouple 
and separate the perturbative equations in an invariant way [15]. In the literature 
there are mathematical studies of the perturbations of type D spacetimes in terms of 
Hertzian and Debye potentials [25-28]. Here we prefer to attack the problem by using 
the standard Teukolsky theory instead. The TME admits separable solutions of the 
form 

xl;'^'\t,x,y,z) = e'''''^e'''^ye'^''-Y{t) , (3.9) 
with Y{t) satisfying the "Master" equation 



dt 



ti^+h)±Y{t) 
dt 



k^,t ^ - 2iskrt 3 + khs 



Y{t) = 0, (3.10) 



where fc^ = ky + k"^. We point out that this ODE can be reduced to a Heun biconfluent 
hypergeometric equation. In fact, introducing the rescaling 

r(t) = t2/3«e* 3 Z(t) (3.11) 
and then changing the variable t as follows 



2 



3k 

Eq. (|3.1U|I takes the canonical form [29,30] 

'^ + (l + a^(3iy- 2i/2) + f (^j - a - 2) v - ^ [S + (1 + a) f3][ Z = , (3.13) 
dv"^ dv \ 2 ' 

with 



e'* v^'^, (3.12) 



V 



9/2 



1/2 

/3 = 0, 7 = 3s, 5 = -^(^-l-j kle^-'\ (3.14) 

Additional perturbative sources can easily be introduced in the problem once expanded 
on this complete basis. Or even more simply, our equation could be easily solved 
for asymptotic values i — > (although in this limit the curvature invariants become 
infinite and perturbation theory clearly has problems) and i — > oo. In analogy with 
black hole physics in which, because of the Peeling Theorem, one has that the fields 
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asymptotically in space are described by radial power laws, here the physical field 
components are described by gauge and tetrad invariant temporal power laws which 
are dependent on s, reminding "a Cosmological Peeling-ofF Property of Gravity" [31]. 
However in the following we will not be interested in this more formal task, but instead 
we will analyze a very simple subset of perturbations, making a link with Section 2. 

We now analyze those special solutions of the TME that represent Kasner models 
perturbatively close to the type D {pi = —1/3, P2 = P3 = 2/3) background one. First 
we write again ps — 2/3 + £ and, taking into account the parameter constraints H2.5|l . 
expanding the expressions (|3.3|) we obtain that the (real) exact Weyl scalars become 

1 4 2 1 

^o^^t-h + o{e^), ^p2~^-t-^ + o{e^), 7/.4 =i 3i"^£ + o(£2) (3 -^5) 

showing again that despite the PNDs are not analytic on the type D point, the 
curvature is. 

Going back to our Master equation H3.10|l . the substitution oi — ky — kz ~ 
in H3.10|l gives the general solution Y = ci+C2t~^^ , which for s — ±2 yields exactly the 
linearized Weyl scalars H3.15|l above. Not surprisingly, to recover a Kasner model from 
the general perturbation equation H3.10|l we have to assume homogeneity, represented 
by vanishing wave- numbers {i.e. infinite wavelengths). 

We point out that for perturbed black holes ipo and ip4 describe gravitational 
radiation at infinity [3]. More in general, using the interpretation provided by 
Szekeres's gravitational compass [9], ^po and ?/'4 are responsible for the transverse 
deformations with respect to the congruences I and n, including in particular the 
possible deformations of gravitational wave type. For the case of our special 
perturbations the wave vector is zero (A: = 0), implying the absence of propagation of 
gravitational signals for perturbations that have to represent Kasner solutions close to 
the LRS type D solution. The absence of gravitational waves in Kasner models is well 
known, and had to be expected here from the vanishing of the magnetic Weyl tensor 
H2.13|l and from the manifest spatial homogeneity of the metric. Although the metric 
tensor is not a gauge invariant object, the spatial independence of (|2.4|) directly reflects 
the invariant spatial Killing symmetries embodied in the spatial tetrad vectors (|2.11() . 
However, the Weyl scalars V'o^^ and ■04^'' that we have derived from the Master equation 
H3.1()|l are gauge and tetrad invariant perturbations of the type D background, so that 
the absence of spatial propagation is a physical information derived here in a novel 
way, using the formalism of the Teukolsky equation for gravitational perturbations. 

Finally, going back to Szekeres's gravitational compass [9] , the electric Weyl tensor 
H2.12|l represents the only direct curvature contribution to the Jacobi (or in particular 
the geodesic deviation) equation, and for any Petrov type / field and any transverse 
frame can be re-expressed as in equation (|2.14|) . It is actually this expression for E(u)ab 
that justifies in general (and not just in a perturbative context) the "transverse frame" 
terminology [6,9]: for a generic tetrad with -01 ^ and/or ■03 7^ there would also be 
longitudinal contributions to (|2.14|) [9], and the magnetic Weyl tensor H{u)ab would 
be non zero, clearly a "tetrad gauge" effect. For type D spacetimes, observers using 
a canonical tetrad (and associated orthonormal one) don't measure any transverse 
contribution. For the general Kasner case, either considering the preferential observer 
and associated tetrad H2.11|l with the corresponding null tetrad (|2.6|1 - H2.7|I . or any 
other null tetrad like (|3.1|l . one has purely time dependent monotonically decaying 
Weyl scalars. Thus, the absence of gravitational waves is once again manifest. On 
the other hand, the transverse contributions 0o and ^4 to E{u)ab in H2.14|l are a good 
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example that the presence of these transverse terms does not necessarily imply the 

presence of gravitational waves. An analogous results can be obtained in studying the 
spacetime of stationary axi-symmetric rotating neutron stars [32]. 

4. Conclusions 

In this paper we have considered quasi-Z) Kasner models that can be parametrized 
as linear perturbations of the special Petrov type D case. The simple calculations in 
Section 2 allow us to conclude that in general a perturbed spacetime will not retain 
the speciality character of the background and therefore will be of general Petrov 
type /. The fact that the indicators of speciality like S or 5norm (which is more 
significant in the perturbative context) may or may not be linear in the perturbations 
only reflects the way they are defined. Thus the fact that S only changes at second 
order docs not necessary imply Petrov speciality for the linearly perturbed spacetimes. 
In fact we have obtained, as a concrete example, four distinct PNDs for the quasi-i* 
Kasner solutions. In other words, in the case of the SI S, we have to consider it 
as a quantity in which the first non vanishing perturbative term has to be retained, 
to be consistent with first order solutions of Einstein equations. This is not too 
surprising, as S is not a dynamical variable, but instead by construction is higher 
order, as the curvature invariants I and J. It is the solution of Einstein equations 
that has to admit a Taylor expansion around the given background [15-17,33], and 
it is in solving Einstein equations at first order that higher order terms have to be 
neglected for consistency. More precisely, the existence of a parametrised family of 
solutions admitting a Taylor expansion around a given background is a prerequisite for 
the application of perturbation methods to Einstein equations [15-17,33]. Once the 
perturbative solution has been found, this does not imply that any quantity of interest 
must be necessarily computed retaining first order terms only. If the first perturbative 
terms in such a quantity are of higher order, the first non vanishing perturbative order 
must be retained. Also, if a Taylor expansion of a spacetime family exists aroimd a 
given background, this does not imply that any quantity computed for the perturbed 
spacetime is a Taylor expansion: as shown by our expansion of the PNDs, they do not 
admit a Taylor series around the type D background case. 

In this article we have also adopted the Teukolsky Master equation formalism, 
which is commonly used in black holes perturbation theory to study gravitational 
radiation. Importing this tool in cosmology we have used it here to demonstrate that 
a quasi-D Kasner spacetime does not contain gravitational waves. This confirms, with 
a novel approach, the well known absence of gravitational radiation in general type I 
Kasner spacetimes, shown here by the vanishing of the magnetic Weyl tensor H{u)ab 
and also resulting by the vanishing in Kasner models of the gravitational superenergy 
flux [4]. 

Finally, we have used this Kasner example to argue that the presence of transverse 
curvature contributions in the Jacobi or geodesic deviation equation (represented 
by the Vo and ^4 terms in the Electric Weyl tensor E{u)ah) does not necessarily 
imply the presence of gravitational radiation. Clarification of the interpretation 
of different transverse frames, to the end of choosing those that only contain the 
gravitational waves degrees of freedom in the transverse part of the curvature at large 
distances from the source, is essential to the problem of wave extraction from numerical 
simulations [34]. 
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